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The present knowledge on charm and leptons is reviewed including the topics discussed 
at this conference and the progress made since the last Quark Matter conference. Special 
emphasis is placed on J/ip production at the SPS which is one of the highlights in the 
field. 

1. Introduction 

Observables related to charm and leptons are of great relevance in heavy ion collisions 
considering the significant information that they have revealed over the last years. The 
low mass dielectron yield measured by the CERES collaboration has provided evidence for 
in-medium effects on the p meson. The J/ip suppression pattern observed by the NA50 
collaboration is believed to be an important signature of the QGP. The intermediate 
mass dimuon excess in A-A reactions at SPS has been first attributed to an increase of 
the open charm cross-section. At the beginning of this year, the first measurement of 
the open charm cross-section in heavy ion collisions has been carried-out at RHIC by 
the PHENIX collaboration. On the other hand, the actual mechanism for charmonium 
production in heavy ion reactions is a subject of intensive current debate. Indeed, the new 
scenarios which have been recently proposed differ totally from the standard approach in 
their predictions. At RHIC they imply a significantly larger J/ ip yield than that resulting 
from the standard suppression scenario. This makes the first J/ip measurements from 
PHENIX of crucial importance. 



2. Low mass dileptons at SPS 

Studies of low mass dielectrons at the SPS continue to attract a great interest since the 
observation of an excess beyond the expected sources in S-U reactions at 200 AGeV [ 1 
and in Pb-Au reactions at 158 AGeV [ 0. The enhancement is located in the invariant 
mass range 0.2 < m < 0.7 GeV/c . It is found to be more pronounced at low p t and 
to increase with centrality. The final analysis of the 40 AGeV Pb-Au data presented at 
this conference [ |3| confirms the preliminary results reported at the last Quark Matter 
conference [ |J: the data taken at 40 AGeV exhibit an even larger enhancement (5.1 ± 
1.3(stat) ± l.O(syst)) than at 158 AGeV (2.9 ± 0.3(stat) ± 0.6(syst)). This is probably 
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due to the maximum baryonic density reached in Pb-Au collisions at the lower beam 
energy. Note that, within errors, the measured dielectron yield is, in the considered 
invariant mass range, similar at both beam energies such that the larger excess observed 
at 40 AGeV actually results from a lower magnitude of the sum of known hadronic decay 
sources. The enhancement can be well accounted for by considering reduction in mass or 
broadening in width of the p meson in the hot and dense medium (see [ |5| for a review 
and [|6| for a recent update). Unfortunately, the present experimental resolution does not 
allow to disentangle the two scenarios (this might be improved with the analysis of the 
data taken in 2000). Another approach, recently proposed, is based on dileptons from 
a simple parameterization of a thermal source [ 0. This approach is very successful in 
describing also the intermediate mass dilepton yield and the direct photons at SPS (see 
below). The derived space-time averaged temperature is 170(145) MeV at 158(40) AGeV. 
It is interesting to note that the temperature at 158 AGeV coincide with the expected 
temperature of the phase transition. Its maximum is located at 210 MeV suggesting that 
the system originates from the deconfmed region [ [7| . 



3. Intermediate mass dileptons at SPS 

In the intermediate mass range (IMR, 1.5 < m < 2.5 GeV/c 2 ), like in the case of the low 
mass region, an excess of dileptons is seen compared to known sources in A-A reactions. 
This has been observed by the HELIOS/3 [ffl and NA50 [ 0] collaborations. An updated 



analysis was presented at the last Quark Matter conference [ [10]. The enhancement is 
found to increase linearly with N part and reaches a factor ~ 2.1 for the most central Pb-Pb 
events. Surprisingly, the kinematical distributions of the dimuon excess are compatible 
with those expected from open charm decay [ ||. These findings have triggered a lot 
of theoretical activities with various interpretations: open charm enhancement [ D 
meson final state rescattering [ [12|, secondary hadron production [ secondary Drell- 
Yan production [ TJ], rj c production [ 115), thermal dilepton production [ M, secondary 



meson-baryon interactions [ jlgl , in- medium Drell-Yan production [ [19) . Some of the 
proposed scenarios fail in reproducing the kinematical distributions of the excess, some 
fail in its centrality dependence and some underpredict its magnitude. The most successful 
description of the data has been achieved by means of thermal dileptons [ [Tj| (see also [ 
T7j). Indeed, with a common thermal source one can simultaneously reproduce the low 



mass dielectron yield from CERES, the IMR dimuons from NA50 and HELIOS/3 as well 
as the direct photons from WA98 [ [7|. However, it should be mentioned that a recent 
investigation of the centrality dependence of the IMR dimuon excess leads to a source 
temperature which is significantly larger than that estimated from central events only [ 



20]] . This could leave room for a partial contribution from one or some of the other effects 
listed above. 

One should finally note that, from the experimental point of view, this kind of analysis 
requires a thorough understanding of the combinatorial background. This is because, in 
contrast to narrow resonances, the present broad signal does not clearly stick-out of the 
background. In particular, special attention should be paid to situations with a production 
asymmetry between positive and negative leptons and a large signal which will lead to an 



incorrect background subtraction and therefore to a biased extracted signal [ [21 
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4. Open charm at RHIC 

The first (indirect) measurement of the open charm cross-section in heavy ion colli- 
sions has been carried-out recently by the PHENIX collaboration via analysis of single 



electron spectra in Au-Au at ^fs = 130 GeV [ 22|. After subtracting all known sources 
from the total distribution the remaining electron spectra are, in the explored p t range, in 
good agreement with electron spectra from open charm decay as predicted using PYTHIA 
(tuned to reproduce existing hadron-nucleon data at lower beam energies and then ex- 
trapolated to Au-Au at RHIC). The derived integrated cross-section is in good agreement 
with PYTHIA and consistent with pQCD NLO calculations. Therefore, in contrast to 
the SPS IMR dimuon data, there is no need for charm enhancement in order to explain 
the RHIC single electron spectra. However, because of the present large systematical 
errors, one cannot firmly exclude in-medium effects on charm production at RHIC. It 
was actually pointed-out at this conference that the small differences at high p t between 
the minimum-bias spectrum and the central spectrum could be explained by considering 



energy loss of charm quarks [ 23 



The analysis was repeated for the 200 GeV data [|24|. Here, the determination of the 
background was done in a completely different way, namely by using a photon converter 
instead of the cocktail technique which was used for the 130 GeV data. Again, the data 
agree well with the expectations from PYTHIA. This holds true for several centrality bins 
such that, here again, not much room is left for large in- medium effects. However, it is 
clear that the safest way to make a definite statement about in-medium effects on open 
charm at RHIC is to compare the Au-Au data to the p-p data. This is hopefully going 
to be done soon [ |24| . 

Another interesting perspective, with the coming high luminosity runs, is the indepen- 
dent measurement of the open charm cross-section from the dielectron continuum above 
the meson. All these measurements could provide precious information on open bottom 
by looking at high p t and high invariant mass. They would be greatly improved by means 
of high resolution vertex detectors owing to the large cr of open charm and open bot- 
tom hadrons. Note finally that non-direct measurements of open charm (bottom) from 
semi-leptonic decays bring limited information only. Direct measurements of open charm 
hadrons remain an important goal of the heavy quark physics program at RHIC. 

5. Charmonium at SPS 

After many years of intensive experimental and theoretical investigations, charmonium 



production at SPS remains one of the hottest topic in the field (see [ |25[ for reviews). 
This is because the observed behaviour of J/ip is believed to exhibit one of the most 
pronounced deviation with respect to pure conventional hadronic picture. 

The so-called "J/ ip anomalous suppression" has been observed by the NA50 collabora- 
tion with the data collected in 1996 and 1998. It manifests itself by a departure from the 
"normal" nuclear absorption when looking at the centrality dependence of the (J/ip)/ (DY) 
ratio in Pb-Pb reactions [ ^6], 27 . The ratio is obtained by means of either the standard 



method or the minimum-bias method. The centrality of the reactions is measured with 
the neutral transverse energy (E t ) and the Zero Degree Calorimeter energy (Ezdc)- A 
~ 20% drop is observed at E t ~ 40 GeV, followed by a decrease at E t ~ 100 GeV (Tab. |). 



4 



The NA50 collaboration interprets this pattern as an evidence for deconfinement. In 
this scenario, the drop results from the melting of XcQ and the decrease from the melting 
of J/ip. This can be checked [ [31] by comparing the energy densities derived from the 
data (Tab. |IJ) with the expected charmonia dissociation temperatures from recent studies 
based on the heavy quark potential from LQCD [ [32]| . The data agree with the predictions 
within ~ 20%. This is a fairly good agreement considering the uncertainties in the way 
the energy density is derived from the data. Note however that these studies assume the 
Xc decay branching ratio into J/ip in heavy ion reactions to be the same as in p-p. 

Table 1 

Location of the drop and the decrease with the transverse energy (E t ), the Zero Degree 
Calorimeter energy (Ezdc)j the path length through matter (L), the impact parameter 
(b), the number of participant nucleons (N part ) and the energy density (e). The values 
are from the data collected in 1996 and 1998 [ |(| [27[| . 





E t (GeV) 


Ezdc (TeV) 


L (fm) 


b (fm) 


N part 


e (GeV/fm 3 ) 


drop 


41 


25 


8 


8 


122 


2.3 


decrease 


100 


9 


9.3 


3 


334 


3.1 



5.1. (J/^)/(DY) in central Pb-Pb 

The interpretation of the decrease as due to the melting of J/ip is not straightforward 
in a dynamical picture. Indeed, QGP models which assume two sharp thresholds corre- 
sponding to the successive meltings of Xc and J /if) overpredict the data in the decrease 
(i.e. beyond E t = 100 GeV) 
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In fact, it was realized that at E t = 100 GeV the 
E t distribution exhibits the typical knee beyond which there is no further sensitivity to 
centrality and fluctuations set in [ [34| . By taking into account these fluctuations, a perfect 
fit of the data is achieved over all the centrality range with either two thresholds or one 
gradual threshold 
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In other words, in the framework of QGP models, the decrease 
seen in the data cannot be non-ambiguously associated to the melting of J/ip. 

Furthermore, it has recently been pointed-out that the decrease may hide another effect 
which has strictly nothing to do with the QGP but which is related to the trigger of the 
NA50 spectrometer 
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An event which fulfills the J/ip trigger conditions should have, 
in average, a smaller E t than a minimum-bias event because a fraction of E t is taken away 
by the J/if>. Therefore, in the minimum-bias analysis, the (J/i/j)/(DY) ratio is biased 
because the J/i[) yield is measured (under trigger conditions) whereas the Drell-Yan yield 
is calculated (without trigger conditions). It has been shown in [ [35| that this effect leads 
to an extra depletion of the (J/V , )/(DY)| min _bias ratio in the E t range of the decrease. 
Consequently, with E t fluctuations and E t loss effects, the comover model provides a very 
good description of the data up to the most central events [ [H| . This also holds true for 
other models based on different scenarios [ BBj, 1371, BB[ and would probably also hold true 



for former models which assume or not the formation of the QGP [ |39j. In this context 
the decrease cannot be considered any longer as a valid point to disprove a model which 
would exhibit a saturation of the (J/^)/(DY) ratio at high E t . Note that the trigger E t 

2 The fraction of Jhj> from Xc decay is ~ 30% in p-p at -y/s = 1.8 TeV [ |§ and - 32% in p-A at 
v / i= 41.6 GeV [M. 
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loss effect can, in principle, be confirmed or disproved by taking Drell-Yan data at large 
E t and/or by measuring at RHIC the J/ if) yield at the energy density corresponding to 
that of the decrease at the SPS. 

5.2. (J/if))/(DY) in peripheral Pb-Pb 

At the time of the last Quark Matter conference, the (J/ip)/(DY) data at low E t were 



subject to different interpretations. Indeed, as it can be seen from Fig. 8 of [ f40|, at low 
E t , some models underpredict the data, some overpredict the data and some go through 
the data. On the other hand, the data points below E t = 30 GeV are located above 
the nuclear absorption curve which means that they show less absorption than in p-A 
and S-U systems. The same feature is observed at large Ezdc [ |27| (i- e - f° r the same 
centrality class). This is difficult to understand in the QGP picture since the suppression 
should be "normal" below the (first) threshold. However, it was believed that, for these 
peripheral events, the data could be contaminated by beam-air interactions. In order to 
clarify the situation new data were taken in 2000 with the target placed in the vacuum 
to minimize off-target interactions. The software was also improved with new methods 
for rejection of pile-up events and identification of the interaction point. In parallel, more 
accurate measurements of p-A data were made in order to estimate with better precision 
the J/if> nuclear absorption cross-section. The updated value is Ojfe = 4.4 ± 0.5 mb [ 
(the former was offo = 6.4 ± 0.8 mb [ 0). The new (J/^)/(DY)| standard data have 
been presented at this conference [ They are consistent with the previously published 
data for E t > 30 GeV. On the contrary, for E t < 30 GeV, they lie significantly below the 
former data and are compatible with the (new) nuclear absorption curve. Consequently, 
models going through the previous data points at low E t should strongly deviate from the 
new data in this E t range. In particular, this is the case for the comover model [ [44|, 55). 



Conclusion: It has been shown at this conference [ |43| that the published peripheral 
(J/ijj)/(DY) data suffer from a systematical bias due to beam-air interactions (a ~ 25% 
effect at E t = 15 GeV). It has been argued in [ ^] that the central (J '/ 'if)) / (DY)| min _ bias 
data might be affected by a systematical bias due to trigger E t loss (a ~ 20% effect at 
Et = 120 GeV). Very recently, it has been claimed that there might be systematical 
inconsistencies between the E t data and the Ezdc data for semi-central events [ 28|| . As 



interpretations are often being drawn at a few percent level, it is extremely important 
that the data are carefully checked for systematical effects. In addition, more stringent 
tests of the models could certainly be achieved i) by studying simultaneously the central- 
ity dependence of yields together with ratios and ii) by extending investigations to all 
available observables. The latter should include, in particular, the transverse momentum 
dependence of the suppression discussed below. 

5.3. Transverse momentum dependence of the suppression 

The p t dependence of the J /if) suppression was early considered to be the golden ob- 
servable to prove the existence of the QGP. Traditionally the data are presented in terms 
of the centrality dependence of the < p\ > of J /if) in order to provide evidence for the 
Pt broadening effect due to initial state parton scattering which is expected in absence of 
QGP. In a deconfined medium, it was first believed that one should observe additional 
p t broadening with increasing centrality [ |45j . This is because, due to the resonance for- 
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mation time and screening effects, only high p t J /if} could escape the medium whereas 
low p t J / if) would melt. Later on, it was claimed [ [46|] that an opposite effect should 
actually happen, namely < pi > flattening or even decreasing. The reason is that high p t 
J I if) are those which travel through the largest amount of matter and therefore are those 
which come from the center of the reaction. Since the local density is the highest in the 
central region these high p t J /if) should be the most QGP suppressed. On the other hand, 
comovers are expected to produce a weak < pi > flattening [ f|7|, [48| . 

The NA50 Pb-Pb data [ || are shown in Fig. 1 (see [ || for the data taken in 2000 
and [ ^TJ for the S-U data at 200 AGeV). They show a smooth increase of < pi > with 
centrality and a saturation when approaching mid-central reactions. It is astounding to 
note that the data do not show any significant sign for onsets, steps or drops and it would 
be remarkable that the two QGP effects mentioned above would exactly cancel each other. 
Even more striking is the outcome of the comparison with the theory, since the only model 
which significantly misses the data is the QGP model (first noticed in [ |52||). This means 
that something is inconsistent when looking simultaneously at {J /if>) / (DY) and at < pi > 
of J /if>. It is therefore of crucial importance to further investigate these data carefully. 
To this respect, it is also important to stress that one would surely learn much more from 
the spectra themselves instead of limiting the study to mean values^. 



CM 



o 

CM 

> 1.9 
<D 

O 



1.8 - 



A 

CM 

CL 

V 1.7 



1.6 
1.5 
1.4 




S. Gavin and R. Vogt [47] 

N. Armesto et al., abs. only [48] 

N. Armesto et al., abs.+comovers [48] 

— D. Kharzeev et al., confined [46] 

■ ■■ ■ D. Kharzeev et al., deconfined [46] 



20 40 60 80 100 120 140 



Figure 1. < T>1 > of J/ip 
as a function of the transverse 
energy in Pb-Pb reactions at 
158 AGeV. The data (dots) [ 
49| are compared to various 



model predictions (curves) with 
references reported on the fig- 
ure. Some of the calculations 
were performed at a beam en- 
ergy of 200 AGeV instead of 
158 AGeV. In order to account 
for this, the < pi > values have 
been rescaled according to [ [52 
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Another interesting aspect of the centrality dependence of the J / if) pt spectra has been 
presented at this conference [ 50|. The J/if> inverse slope parameters for p-A, S-U and non- 
central Pb-Pb reactions fall on a common straight line when plotted as a function of the 
energy density. However, for the most-central Pb-Pb events one observes a flattening of 
the apparent temperatures. This could indicate a change in the expansion of the system. 

5.4. Non-direct charmonium production 

Charmonium production from statistical hadronization has been extensively discussed 
in a parallel session at this conference. This production mechanism is based on a twofold 

3 As illustrated in [ |53|, the fact that a model gets the mean of a distribution right does not necessarily 
mean that the full distribution is understood. 
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observation [ [54 . First, because of cross-section considerations, the treatment of charmed 
hadrons in a standard thermal model cannot be easily justified. Secondly, in central Pb-Pb 
reactions at SPS the measured if)'/ ( J /if)) ratio is ~ 4% as predicted by the thermal model. 
This led to a new scenario for charmonium production [ |54j] : i) all cc pairs are produced in 
the early stage via hard scattering; ii) all charm hadrons are produced statistically at the 
hadronization. The underlying picture is that the charm hadrons either do not form in 
the early stages or are fully suppressed in the QGP. This scenario can be implemented in 
the framework of a thermal model by using the formalism described in [ |)4j |55|, [56] where 
one includes a charm enhancement factor in order to match the multiplicity of directly 
produced cc pairs. Since the open charm cross-section has not yet been measured at SPS, 
it is taken from the predictions of pQCD. 



This approach has been studied in detail in the two last years [ 54], [56| [571 |36L |60| 



pT| , p7[ . As shown in [ |5B[ the measured J /if) yield versus N part in Pb-Pb is qualitatively 
consistent with the centrality dependence of the nuclear overlap function. In order to 
reproduce the yields, one needs to enhance the charm cross-section in the model with 
respect to the predictions from pQCD. This extra charm enhancement amounts to a factor 
~ 2. This is surprisingly similar to the open charm-like enhancement needed in order to 
account for the excess of the IMR dimuon yield observed in central Pb-Pb events at SPS [ 
|Kj. Another realization of the statistical hadronization model also provides a very good 
fit of the E t dependence of the (J/if>)/ (DY) ratio, albeit with a larger enhancement factor [ 



57, >A 



An additional experimental indication for statistical hadronization of charmonium at 
SPS has been put forward in [ |58]] (see [ for the model predictions at RHIC): The m t 
spectra of Q ± from WA97 and J /if) and if)' from NA50 can be simultaneously described in 
a hydrodynamical picture of hadronizing QGP with T = T c = 170 MeV and v = 0.2 (v 
is the average transverse flow velocity). This supports the idea of charmonium formation 
and kinetic freeze-out at the hadronization. 

The statistical hadronization model can also be used to predict charmonium yield at 
RHIC [ [5IJ. At RHIC, depending on the beam energy, the total multiplicity of (directly 
produced) cc pairs is expected to be smaller than one in peripheral reactions and larger 
than one in central reactions. Because the statistical hadronization process is strongly 
correlated to the multiplicity of the cc pairs, one expects the final (J/ip)/(N c5 ) ratio to 
increase or decrease versus centrality depending on the beam energy [ ^UJ . 

The J / if) statistical hadronization process has been recently incorporated in a modified 
approach [ ^lj, [37|. This approach combines i) directly produced J /if) followed by nuclear 
absorption and QGP/HG dissociation and ii) statistically produced J /if) as in pure sta- 
tistical models but without the extra charm enhancement factor. This leads to a good 
agreement with the NA50 (J/if))/(DY) data. The agreement is particularly good at the 
level of the first drop which arises, in this model, from the superposition of direct and 
statistically produced J /if). Another interesting aspect from this model is the excitation 
function of the total J / if) yield which is expected to be dominated by direct J / if) at SPS 
and by statistical J /if) at RHIC, with a minimum at y^i ~ 50 GeV. It is important to 
note that whereas pure statistical models reproduce the measured if)' /{J /if)) ratio, this 
approach significantly overpredicts it. This leaves room for additional effects [ |6T| . 

Non-direct charmonium states can be produced not only at the hadronization but also 
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in the QGP by means of the so-called kinetic recombination process [ ^2] . This mechanism 
foresees even larger J/ip yields at RHIC than the statistical hadronization process. Note 
that all models will face stringent tests with the upcoming data from NA60 at SPS and 
PHENIX at RHIC. 

To conclude on this topic, it is interesting to extrapolate from charmonium production 
at RHIC to bottomonium production at LHC. Indeed, i) the key parameter for statistical 
production of charmonium is the number of cc pairs and ii) the expected multiplicity of bb 
pairs at LHC is roughly equal to the expected multiplicity of cc pairs at RHIC. Therefore, 
if statistical production of charmonium is confirmed at RHIC, one could expect the same 
mechanism to occur at LHC for bottomonium states. 



6. Charmonium at RHIC 

The first charmonium measurements from the PHENIX collaboration are certainly 
among the most exciting results presented at this conference [ ^3[. The J/ip signal has 
been extracted in p-p at \/s = 200 GeV both in the dielectron channel and in the dimuon 
channel. Although the present statistics is very limited (~ 24(36) counts in the dielec- 
tron(dimuon) channel), the signal is clearly seen. Thanks to the capability of the PHENIX 
detector to measure simultaneously the two channels, the first rapidity distribution of J / ip 
in p-p at y/s = 200 GeV has been established. The J/ip signal has also been extracted 
in Au-Au reactions in the dielectron channel. The statistics, even more limited than in 
p-p, does not allow to draw at present solid conclusions on a possible J/ip suppression 
or enhancement. However, large suppresion factors as well as large enhancement factors 
seem to be not likely. This will be improved soon since only half of the collected data 
from the 200 GeV run has been analysed so far. In the next two years the statistics will 
be enhanced by a factor ~ 75 for Au-Au and ~ 100 for p-p. Also of great interest will be 
the coming charmonium measurements in d-A. 



7. Next steps 

Apart from the forthcoming measurements at RHIC and the analysis of remaining SPS 
data, new generation experiments are presently being completed or designed. 



The HADES experiment [|6"3| at GSI will soon shed light on the DLS low mass dielectron 
puzzle [j65|. The NA60 experiment [^6| at SPS will continue in line with the investigations 
made by NA50 and, in particular, measure the open charm cross-section and the ip' . The 



ALICE experiment [ p7| will bring essential information regarding charm and leptons 



at LHC. Looking even forward into the future, the CBM experiment [ |68|, at the GSI 
future facility, will allow a new beam energy domain to be explored with simultaneous 
measurements of low mass dielectrons, charmonia and open charm. 
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